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We study soliton pulse compression in materials with cascaded quadratic nonlinearities and show that 
the group- velocity mismatch creates two different temporally nonlocal regimes. They correspond to 
what is known as the stationary and nonstationary regimes. The theory accurately predicts the transition to 
the stationary regime, where highly efficient pulse compression is possible. © 2008 Optical Society of America 
OCIS codes: 320.5520, 320.7110, 190.5530, 190.2620, 190.4400 



Efficient soliton pulse compression is possible using 
second-harmonic generation (SHG) in the limit of large 
phase mismatch, because a Kerr-like nonlinear phase 
shift is induced on the fundamental wave (FW). Large 
negative phase shifts can be created, since the phase mis- 
match determines the sign and magnitude of the effective 
cubic nonlinearity. This induced self-defocusing nonlin- 
carity thus creates a negative linear chirp through an ef- 
fective self-phase modulation (SPM) term, and the pulse 
can therefore be compressed with normal dispersion. 
Beam filamentation and other problems normally en- 
countered due to self-focusing in cubic media are there- 
fore avoided. This self-defocusing soliton compressor can 
create high-energy few-cycle fs pulses in bulk materials 
with no power limit [1-4]. However, the group- velocity 
mismatch (GVM) between the FW and second harmonic 
(SH) limits the pulse quality and compression ratio. Es- 
pecially very short input pulses (< 100 fs) give asymmet- 
ric compressed pulses and pulse splitting occurs [4,5]. In 
this case, the system is in the nonstationary regime, and 
conversely when GVM effects can be neglected it is in the 
stationary regime [3-5]. Until now, the stationary regime 
was argued to be when the characteristic GVM length is 
4 times longer than the SHG coherence length [1] , while a 
more accurate perturbative description showed that the 
FW has a GVM-induced Raman-like term [4,5], which 
must be small for the system to be in the stationary 
regime [4]. However, no precise definition of the transi- 
tion between the regimes exists. 

On the other hand, the concept of nonlocality provides 
accurate predictions of quadratic spatial solitons [6,7], 
and many other physical systems (see [8] for a review). 
Here we introduce the concept of nonlocality to the tem- 
poral regime and soliton pulse compression in quadratic 
nonlinear materials. As we shall show, GVM, the phase 
mismatch, and the SH group-velocity dispersion (GVD) 
all play a key role in defining the nonlocal behavior of 
the system. Two different nonlocal response functions 
appear naturally, one with a localized amplitude - rep- 



resenting the stationary regime - and one with a purely 
oscillatory amplitude - representing the nonstationary 
regime. In the presence of GVM they are asymmetric 
and thus give rise to a Raman effect on the compressed 
pulse. 

In the theoretical analysis we may neglect diffraction, 
higher-order dispersion, cubic Raman terms, and self- 
steepening to get the SHG propagation equations for the 
FW (wi) and SH (uj 2 = 2wi) fields E h2 (z,t) [3,9]: 

(id, - ^dt^E, + Kl E$E 2 e iAk * (1) 
+ Pl E 1 ([E^ + r,\E 2 \ 2 ) =0, 
(id, - id 12 d t - \kfd tt )E 2 + K 2 E 2 e~* Akz 

+ p 2 E 2 (\E 2 \ 2 + r 1 \E 1 \ 2 ) =0, (2) 

where r\ = 2/3 for Type I SHG [3]. Kj = uj\d e s/ arij, c? ff 
is the effective quadratic nonlinearity, pj = idjUKen j/c, 
and riKcrr.j = 3Re(%' 3 ) )/8n,' is the cubic (Kerr) non- 
linear refractive index. The phase mismatch is AA; = 
k 2 — 2fci, with kj = UjUjj / c. rij is the refractive index, and 
kj = d n kj j l du n \ u>=UJj accounts for dispersion. The time 
coordinate moves with the FW group velocity v s ,i = 
l/kr[ \ giving the GVM term d\ 2 = v~{ — v~^. In dimen- 
sionless form, r = t/Tb, where To is the FW input pulse 
duration, £ = z/Ld,i, where Ld,i = To/\k^\ is the 
FW dispersion length, and U\ = Ei/Eo, U 2 = E 2 j\fr\£o, 
where Bo = Ei(0,0) and n = ni/n 2 . Equations (TTJ) and 
© become [3] 

(id^ - D 1 d rr )U 1 + ^XpNs W UlU 2 e lAK 

+ Nl m U 1 (\U 1 \ 2 + V n\U 2 \ 2 )=0, (3) 

(id e - tSd T - D 2 d TT )U 2 + ^A0N SiiG U 2 e- lA ^ 

+ 2n 2 Nl m U 2 (\U 2 \ 2 + nh- 1 ^ 2 ) = 0, (4) 

with S = di a To/|*i 2) |, Dj = kf ) /2\k ( i 2) \, and A/3 = 
A/cLd,i- The scaling conveniently gives the SHG soliton 
number [3,4], V| HG = Lv.i£qUj 2 d 2 s / (c 2 nin 2 Ak), and 
the Kerr soliton number Nfi_ elI = ^D,i^Kerr,i^o^i/c 
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In the cascading limit Aj3 ^> 1 the nonlocal approach 
takes C/ 2 (£,r) = </> 2 (r) exp(— «A/3£), keeping its time de- 
pendence but neglecting the dependence on £ of </> 2 . 
To do this the coherence length L co ^ = 7r/|Afc| must 
be the shortest characteristic length scale in the sys- 
tem, which is true in all cascaded compression exper- 
iments. Discarding the Kerr terms in Eq. (j4|) because 
N 2 arv U 2 <C \/A(3Nshg, we get an ordinary differential 



Kerr 

equation D 2 (f> 2 + it 
formal solution 



Mr) = - 



A^SHG 



•A/302 = v^TVshg^i, with the 



&t'R ± {t')UI{t-t'). (5) 



According to the sign of the parameter s\ 
sgn (A/3/D 2 - 5 2 /4D%), R + or i?_ must be used 
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sin(|r|/n), 



(6) 
(7) 



where s 2 = sgn(£> 2 /<$) and r^ o dTi?+(r) = 1. Both 
response functions have an asymmetric imaginary part 
due to GVM, causing a Raman effect. Moreover, R + is 
localized in amplitude (corresponding to the stationary 
regime) , while R- is purely oscillating in amplitude (cor- 
responding to the nonstationary regime). The nonlocal 
dimensionless temporal scales are 



n = \A{3/D 2 -t 2 



1-1/2 



t 2 = 2\D 2 /6\. 

.(2) 



(8) 



On dimensional form ti = t x T = \2Ak/k K 2 ' - t 2 2 \-^ 2 , 

(2) 

t 2 = t 2 T = \k 2 /di2\, and TZ± = R±/T , which all are 
independent of Tq. The transition between the stationary 
and nonstationary regimes occurs when si changes sign, 
which in dimensional units implies that 

d\ 2 < 2Akk^ 2 \ stationary regime, R + , (9) 

(10) 



9 (2) 
d 12 > 2Akk 2 , nonstationary regime, R- 



The transition is independent of To but depends on 
the GVM, the SH GVD, and the phase mismatch. This 
central result is important wherever cascaded quadratic 
nonlinear phase shifts are used. It should be compared 
with the qualitative arguments of [1], where the station- 
ary regime was Ak < 47r|di2 |/Tb. 

Now, inserting Eq. ([5]) into ([3]) gives a dimensionless 
generalized nonlinear Schrodinger equation 

^- J D 1 9 TT ](7 1 +iV2 crr [/ 1 |t/ 1 | 2 

/OO 
dT'R ± (r')U*(T-T') = 0, (II) 
-00 

In the weakly nonlocal limit of the stationary regime 
Ti,2 "C 1, U 2 is slow compared with R + , so 

[idt - D 1 d TT ]U 1 - (iV| HG - Nl^WUxf 

= is 2 N 2 ^ G r R ^ G \U 1 \ 2 d T U u (12) 
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Fig. 1. (a) Compression window Eq. (fT4|) versus Ai in a 
BBO. Also Afc sr = 47r|di 2 |/T of [1] for a 100 fs pulse is 
shown, (b) ii j2 and Tr^hg versus Ak for Ai = 1064 nm. 



where f dtR + (t) = 1 was used. The dimensionless 
characteristic time of the nonlocal Raman response is 



T, 
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(13) 



Exactly this result has been derived before by pertur- 
bative methods [3-5]; that method therefore amounts 
to being in the weakly nonlocal limit of the stationary 
regime. Clearly, both Eq. (fTT|) and (fT2|) have terms rem- 
iniscent of the nonlinear Schrodinger equation with a 
purely cubic nonlinearity [10]; the GVM-induced nonlo- 
cality is similar to the Raman terms from a delayed cubic 
response. 

From Eqs. (J9j) and (fTTj) clean soliton compression re- 
quires, first, being in the stationary regime, i.e., Ak > 



Ak s - 



d\ 2 /2k 2 . Second, soliton compression requires 

\J ^shg — ^Korr > 1 P]- This can ^ e expressed 
as Ak < Ak c = Afc c , max /(1 + N^ 2 rl ), where Ak c , max = 
^id 2 s /nKorr,icnin 2 . To remove the dependence on the 
FW input pulse intensity and duration it is convenient 
to require Afc < Afc Cjlnax , which is necessary to have 
iV e fj > 1. Thus, we obtain a compression window: 



Ak. r < Ak < Ak r 



(14) 



In Fig. [TJa) we show the compression window for a 
/3-barium-borate (BBO) crystal (see [3] for BBO mate- 
rial parameter details). Notice that the window closes for 
Ai < 750 nm. Opening the compression window here re- 
quires a material with a stronger quadratic nonlinearity 
or alternatively a strong dispersion control, as offered by 
photonic crystal fibers [11]. At Ai = 800 nm the window 
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Fig. 2. Numerical simulations of soliton compression of 
a 100 fs scch pulse in BBO for Ai = 1064 nm and 
-^eff = 5.3. The normalized FW intensity at z opt is shown 
for (a) Afc = 50 mm" 1 (stationary regime) and (b) 
Afc = 30 mm -1 (nonstationary regime). The response 
functions R±/T and JwnlT) induced by cascaded SPM 
for (a) and (b) are shown in (c) and (e) and (d) and (f), 
respectively. 



is narrow. In fact, the compression experiments done at 
800 nm [1,2] were both in the nonstationary regime, and 
were both unable to observe compression to few-cycle 
pulses. Choosing Ai = 1064 nm, Fig. QJb) shows the 
nonlocal time scales versus Afc. At the transition to the 
nonstationary regime, t\ diverges while £2 remains low. 

Let us illustrate the two regimes. Figures [2ja) and 
[5Jb) show two numerical simulations of soliton compres- 
sion of a 100 fs pulse at Ai = 1064 nm in a BBO (the 
full coupled equations of [3] are solved). In Fig. [2£a) 
Afc = 50 mm -1 ; our nonlocal theory predicts Afc sr = 
36.0 mm , so this is in the stationary regime. Indeed, 
a symmetric compressed 6 fs pulse is observed. Instead, 
changing to the nonstationary regime, Afc = 30 mm -1 
in Fig.[2fb), the pulse at z = z opt (optimal compression 
point) is very asymmetric and strong pulse splitting oc- 
curs. Note, the definition in [1] predicts this simulation 



to be in the stationary regime. The nonlocal response 
functions arc shown in Figs. [2jc) and EJJd); while the 
nonlocal time scales are clearly quite similar for both 
examples, the different shapes of the response functions 
imply a very different impact on the pulse dynamics. 
This can be understood by using Eq. (fTTj) to calculate 
the chirp 5wnl induced by SPM from the cascaded SHG 
process of a sech input pulse (calculations as in [10], 
Chap. 4). Figures (2^0) and [21(f) show <5wnl for z = Lshg 
(where -ZV| HG = £d.i/^shg PD- I n the stationary case, 
Fig. [He), SPM induces a linear negative chirp on the 
central part of the pulse because i?+ is localized, except 
for very short pulses, To ~ t\, where R + becomes non- 
local. In the nonstationary case, Fig. [2ff), SPM induces 
a linear chirp only for long pulses > 30 fs, while shorter 
pulses also have strong chirp induced in the wings be- 
cause of the oscillatory character of i?„ ; this explains the 
trailing pulse train in Fig. [21(b). Thus, the nonstation- 
ary response can be nonlocal even for To 3> ii.2- Note, 
in Fig. [2^b) very short pulses (To ~ £1,2) are positively 
chirped in the central part (equivalent to a chirp induced 
by a self-focusing nonlinearity), making few-cycle com- 
pressed pulses impossible in the nonstationary case. 

To conclude we showed that GVM induces asymmet- 
ric nonlocal Raman responses that accurately explain 
the stationary and nonstationary regimes in cascaded 
quadratic soliton compressors. The nature of the re- 
sponse functions and their degree of nonlocality relative 
to the input pulse length is vital for the resulting com- 
pression. The theory offers new insight into the physics 
of this soliton compressor. 
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